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Разработка программного комплекса 
для оптимизации систем 


В статье рассматривается программный комплекс для оптимизации систем. Комплекс использует 
объектно-ориентированный подход при реализации программ, имеет удобный интерфейс, не требует 
от пользователя знания особенностей языка программирования или деталей схемы организации 
процесса вычислений в комплексе. 


Введение 


Одной из особенностей применения задачи оптимизации к задаче идентификации 
параметров динамической модели, представляющей собой систему дифференциаль- 
ных уравнений, по критерию минимизации несогласованности переходных является 
тот факт, что в таком случае вычисления значения /(х) требуют значительных вы- 


числительных затрат, что затрудняет нахождение на А-й итерации /(х») и делает почти 
невозможным вычисления значения производных /’(х,). Это накладывает определен- 


ные ограничения на класс допустимых алгоритмов оптимизации при решении задачи 
идентификации параметров модели. Данное замечание касается и реальных задач опти- 
мизации технологических режимов функционирования объектов, где в большинстве 
случаев трудно в явном аналитическом виде представить критерий оптимизации. 

Эти соображения стали основой при выборе множества алгоритмов оптимизации, 
которые представлены в разработанном комплексе программ. Каждый из существую- 
щих алгоритмов многомерной оптимизации имеет свои преимущества и недостатки и 
может более эффективно работать или в окрестности точки минимума, или на началь- 
ных этапах оптимизации, когда приближение находится далеко от точки минимума. 

Целью данной работы является разработка комплекса алгоритмов оптими- 
зации для последовательного подключения различных алгоритмов на тех этапах 
процесса оптимизации, где каждый из них будет наиболее эффективным. Существенной 
помощью в выборе эффективного метода в случае конкретной функции /(х) может 
стать предусмотренная в комплексе программ возможность графического представ- 
ления линий уровня исследуемой функции /(х). 


Постановка задачи минимизации 


Разработанная подсистема позволяет определять точку минимума заданной функ- 
ции многих переменных /(х),х = (х1,...,х„) для задачи без ограничений 


аго тт Г (х), 
хЕЁ,„ 


или для задачи с ограничениями 


аге т Л (х), В: а, 5х, «В=Ьп}. 


Ограничения на параметры диктуются самой постановкой задачи, исходя из 
физических соображений. 
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Общая схема реализации алгоритмов минимизации 


В вычислительных схемах большинства алгоритмов минимизации можно вы- 
делить ряд общих последовательных этапов, а именно: 
1) выбор начального приближения хо (нулевая итерация К = 0) и вычисления /(х,); 


2) определение направления спуска р, для минимизации /(х) (способ решения 


этой подзадачи определяющий, для каждого конкретного алгоритма характерна своя 
стратегия при выборе направления спуска); 


3) определение точки минимума х“° на луче: 


шш 
х=х,+а-й,рь: т Ло, +а-й,р,) =, +атйр,); 
4) переопределение полученной точки минимума на луче х,/=х,+ай,рув 


качестве исходной для следующей итерации А =А +1; 

5) условный переход на пункт 2, если критерий завершения (выхода из алго- 
ритма) не выполняется. Однократное выполнение этапов 2 — 5 будем называть одной 
итерацией процесса минимизации. На каждом из этапов нужно фактически решить 
некоторую самостоятельную подзадачу. Для решения подзадач 1 и 3 в комплексе 
предусмотрено несколько альтернативных методов, которые реализованы в виде само- 
стоятельных программных модулей. 

Выбор начальной точки х, осуществляется по такой схеме: если начальная точка 


х, пользователем не задана, то автоматически выбирается и случайная точка из множе- 
ства / =п+1 равномерно распределенных точек &, ‚ К =1,/ в некоторой заданной области 
О, - п-мерном параллелепипеде: /), = |х: А ьи} ‚ где достигается ит Е 


На этапе 5 в понятие «критерий выхода» укладывается проверка выполнения 
двух условий: 

1) алгоритм «исчерпал себя», что отразилось на выполнении так называемого 
«внутреннего» критерия (этот критерий свой для каждого алгоритма), тогда выход из 
программы минимизации происходит с кодом возврата Ку? = 1; 

2) были выполнены заданное число К/Т итераций, максимум на которые был 
запущен данный метод (точнее модуль его реализует), однако внутренний критерий 
не исполнился, тогда возвращение из модуля с кодом Аи = 2. Число итераций К/Т 
задается пользователем или по умолчанию выбирается управляющей программой. 
При возвращении в управляющую программу минимизации с Ки’ = 2 осуществляется 
анализ скорости сходимости алгоритма минимизации. Для этого используется соб- 
ранная в процессе выполнения последних М/Т итераций информация о значении 


1 р } ‚ Е =1, МТ . Эта информация собирается и накапливается в глобальной струк- 


туре 5ТАТ, которая подключена к программным модулям минимизации. Потом идет 
проверка условия К, > ГА |) К, >ГК ‚, где 


ео 
Х МПТ — — А К-+1 2 

1 МТГ-1 ь 
Ю = х.)— Р(х, р 
Е >И Ла) 


Если «нет», то скорость сходимости данного алгоритма неудовлетворительная 
и управляющая программа на следующие К/Т итераций автоматически подключит 
следующий алгоритм минимизации, который предусмотрен в цепочке алгоритмов, 
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которыми планируется решать данную задачу минимизации. Эта цепочка или «по 
умолчанию» назначается автоматически управляющей программой, или может быть 
задана перед началом процесса решения задачи пользователем. Благодаря этому в 
случае неудовлетворительной сходимости пользователь может управлять алгоритмами, 
изменяя значения их входных параметров или конфигурацию модулей, включенных 
в цепочку, и описывать конфигурацию алгоритма. 

Если это условие выполняется, то алгоритм с данной конфигурацией модулей и 
значениями их входных параметров запускается повторно без изменений по следую- 
щим К/Т итерациям. Этот процесс выбора алгоритма на следующее число итераций 
осуществляется до тех пор, пока с алгоритма не произойдет выход по его внутрен- 
нему критерию или заданная цепочка алгоритмов не исчерпает своих возможностей 
по успешной минимизации функции. 

Значения МИГ, ГА, ГК и КИТ задаются пользователем или автоматически опреде- 
ляются «по умолчанию» управляющей программой. Конфигурация алгоритма решения 
задачи, задания значений входных параметров составляющих модулей минимизации, а 
также описание цепочки решающих алгоритмов дается в специальной структуре ТАЁС. 


Множество алгоритмов минимизации 


Для нахождения решения вышепоставленных задач минимизации предусмот- 
рена возможность использования любого допустимого для данной задачи метода, 
реализованного в комплексе программ набора методов по нахождению минимума 
многомерной и одномерной функции. Пользователь имеет возможность по своему 
усмотрению выбирать любой из следующих методов многомерной минимизации: 

1 МСЛН - метод случайных направлений; 

МГРД -— метод градиентного спуска; 

МСМП - симплекс метод; 

МСГР -— метод сопряженных градиентов; 

МПАУ - метод Пауэла; 

МСЛМ - модифицированный метод случайных направлений; 
МОВР - градиентный метод овражного типа; 

МГРС - градиентный метод с растяжением пространства; 
МПКС - метод покоординатного спуска; 

10 МФПД - метод Флетчера — Пауэла — Давидона; 

1] МГРО - метод градиентного спуска с ограничениями; 

12 МСЛО — метод случайных направлений с ограничениями. 

Для реализации одномерного поиска точки минимума в заданном направлении 
комбинировать его с одним из методов одномерной минимизации: 

1 ПАРРТ - метод парабол по равномерным точкам; 

2 ФПАРБ - метод парабол с золотым сечением; 

3 ФИБОН - метод Фибоначчи. 

При решении задачи минимизации целесообразно комбинировать комплексное 
использование алгоритмов на разных этапах алгоритма минимизации в зависимости 
от свойств конкретного алгоритма. С этой точки зрения представленные в подсистеме 
алгоритмы можно классифицировать следующим образом. 


На начальном этапе поиска, когда точка начального приближения Хх, находится 


оон мифьшь 


далеко от точки минимума х„,„, , Целесообразно применять методы линейной аппрокси- 


мации, которые требуют минимальных вычислительных затрат, как, например, градиент- 
ный метод, метод случайных направлений, модифицированный метод случайных 
направлений. Все эти методы имеют линейную скорость сходимости и за сравнительно 
небольшое количество итераций позволяют приблизиться в окрестность точки миниму- 
ма или спуститься на дно «оврага», где скорость сходимости этих методов замедляется. 
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При структуре линий уровня /(х) = с0оп5ё типа «желоба» (вогнутость линий 


будет в направлении малочувствительных компонент) сходимость всех перечисленных 
методов будет плохой. В таком случае возможно «исправление» поверхности путем 
изменения масштабов переменных при применении градиентного метода с растяже- 
нием пространства или применения градиентного метода «овражного» типа. 

При приближении к точке минимума с помощью линейных методов происхо- 
дит замедление сходимости. И поэтому здесь целесообразно переходить на методы 
минимизации второго порядка, использующие квадратичную аппроксимацию /(х). 


Из числа методов такого типа в комплексе представлены: метод Пауэла и метод Флет- 
чера — Пауэла — Давидона. Эти методы минимизируют квадратичную функцию; на 
строго выпуклых дважды непрерывных дифференцируемых функциях имеют квадра- 
тичную сходимость; без предположения строгой выпуклости методы имеют надлиней- 
ную сходимость. Можно также использовать метод сопряженных градиентов, который 
хоть и является методом первого порядка, но имеет более высокую сходимость по срав- 
нению с методами градиентного типа и для квадратичной функции, как и методы Флет- 
чера и Флетчера — Пауэла — Давидона, совпадает не более чем по и итераций. 

Каждый из этих методов реализован в виде отдельного модуля, который может 
быть использован как стандартная подпрограмма из библиотеки модулей и может 
быть подключен к задаче минимизации через посредничество предусмотренного в 
комплексе интерфейса пользователя, что позволяет при решении конкретной задачи 
свести к минимуму этап программирования. 


Организация модулей для задачи оптимизации 


После старта управляющей программы 5"йип() пользователю предоставляется 
возможность с помощью диалоговых процедур задать функцию, которую он хочет 
минимизировать и определить (при желании) алгоритм минимизации. Модуль 5йии() 
осуществляет подготовку к началу процесса решения задачи и инициализирует струк- 
туры: аа — с информацией о поставленной задаче; пой — с описанием модели, в 
рамках которой решается задача; {а/е — с описанием алгоритма решения задачи и 
значениями входных параметров модулей минимизации; ру"! — с описанием режимов 
документирования задачи. Описание этих структур можно найти в файле 16[<[5.срр. 
Далее в модуле яйии() инициализируется заданная функция минимизации мй_тй() 
и происходит задание значений параметров «по умолчанию» и конфигурации 
алгоритма решения задачи. 

За вызов необходимого метода многомерной оптимизации отвечает модуль 
тпитт.срр, методы одномерной минимизации вызываются модулем тотт.срр, 
модуль /[отит.срр отвечает за локализацию точки минимума в заданном направлении. 
Модуль тотт.срр вызывается модулем титт.срр в том случае, если выход из 
подпрограммы минимизации произошел не из-за нахождения точки минимума с задан- 
ной точностью, а через заданное число итераций, т.е. если имеется плохая сходи- 
мость алгоритма. Об этом случае свидетельствует значение параметра Ки = 2. 


Описание процесса разработки модулей 
своей задачи минимизации 


Он включает следующие этапы: 

Создать модуль целевой функции (на основе шаблона функции /()): — вместо Гдать 
имя модуля целевой функции у = Хх). 

В операторе у =... записать формулу для вычисления /(х). Создать модуль опи- 
сания и инициализации задачи (на основе шаблона функции мй{ ()): 
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—ви=... указать размерность вектора х, а в описании сйаг * пате _х[] — наиме- 
нование компонент (для документирования результатов); 

— если задана начальная точка х0, то в описании 4оиЫМе х0/] = () задается ее 
числовое значение и флажок 1.1х0 = 1, иначе под х0 только резервируется память 
аоиЫе х0[п] и 1х0 = 0; 

— если заданная область начальной точки х0, то включаем описание массивов 
аоиЫе 4а0/] = (), 4ои Ме аАБО/] = (), где задаем числовые значения области поиска на- 
чальной точки: 4а0 <х0 <4ЬО, и значение флажком (1.140 = 1; 

— если заданная область ограничений на вектор х, то включаем описание мас- 
сивов доиЫе аа/] = (), 4доиЫе АБ/] = (), в которых задаем числовые значения области 
поиска точки минимум: 4а <х <4Ь, и значение флажком (1.14 = 1. 


/ / «МОДУЛЬ ЦЕЛЕВОЙ ФУНКЦИИЕ (Х)=...», 
дочЫе # (доцЫе * х)(4оцЫе у; у=...; 


теб у; уо1А шит (ТММ &); 
) ТМПХ 6 
/ * МОДУЛЬ инициализации ЗАДАЧИ * / эвайс 4опЫе х0 [] = (); “айс допЫе 4а0 [] = (); 
ТМИХ ( зас 4опЫе 40 [| = (; майс дочЫе аа [|= 0; 
шеи; 16 1х0 = 0; п 140 = 0; шё1а = 0; айс доцЫе 4 [] = (); 
дочЫе * х0; айс саг * пате х []|={«Х [Ш«»Х [2]; 
дочЫе * 4а0, 960; дочЫе * да, 46; фи =; 1х0 =Я.190 =; 414 =; 
сраг * паше []; {-> пате = пате х; 
шит (& 0); 
ус! шИЕ (у01а) теги; 
) 
Заключение 


Разработанный комплекс предназначен для оптимизации систем и благодаря ис- 
пользованию объектно-ориентированного подхода при реализации программ имеет 
удобный интерфейс, не требует от пользователя знания особенностей языка программи- 
рования или деталей схемы организации процесса вычислений в комплексе. 
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В.М. Петрович 

Розробка програмного комплексу для оптимвВацй систем 

У стати розглядаеться програмний комплекс для оптимвацй систем. Комплекс використовуе об’ектно- 
ор1ентований шджд при реаллзацй програм, мае зручний 1нтерфейс, не вимагае в1д користувача знан- 
ня особливостей мови програмування або деталей схеми органзаци процесу обчислень в комплекс1. 


у.М№. Реагомсй 

Оеу@ортеп: о? ЗоЁйумаге ог Орйпи7айоп Зузет 

ТБе агасе 415си$5е$ а зоЙ\аге расКаее Юг орйпияше зузепл$. ТНе сотшрех и5ез Ше офесномеще4 арргоасВ 
40 ргоэтатлилие, Ваз а сопуешепе пщегЁасе фай 4оез пой гедите апу Кпо\Ледее оЁ Фе ргозтатие |апопазе ог 
Фе деа15 ое зсВете$ оЁ сопршаНоп ш Фе сотр/ех. 


Статья поступила в редакцию 08.07.2010. 


70 «Искусственный интеллект» 472010 


